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relativistic  axial  gyrofrequency  is  comparable  to  the  Free  Electron  Laser  Doppler  Upshift, 
in  addition,  there  is  a  range  of  parameters  for  which  the  ponderomotive  potential  acts  to 
destabilize  the  electron  beam.  In  this  regime,  we  find  both  unstable  electrostatic  beam 
modes  and  largely  electromagnetic  modes  with  broad  bandwidths. 
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COLLECTIVE  EFFECTS  ON  THE  OPERATION  OF 
FREE  ELECTRON  LASERS  WITH  AN 
AXIAL  GUIDE  FIELD 


I.  INTRODUCTION 

The  physical  process  which  gives  rise  to  wave  amplification  in  free  electron  lasers  stems  from  the 
interaction  of  a  relativistic  electron  beam  with  a  spatially  periodic  magnetic  field  (i.e.,  the  wiggler  or 
pump  field)  applied  largely  transverse  to  the  direction  of  bulk  electron  motion.  The  effect  of  the 
wiggler  field  is  to  provide  a  coupling  between  the  electron  beam  and  electromagnetic  radiation  fields 
which  results  in  a  ponderomotive  force  along  the  axis  of  the  beam.  The  form  of  the  interaction  can  be 

classified  in  a  variety  of  ways  depending  upon  such  parameters  as  the  magnitude  of  the  electron 

* 

current,  the  strength  qf  the  pump  field,  the  bulk  energy  and  energy  spread  of  the  beam,  and  the  length 
of  the  interaction  region.  For  example,  on  the  one  hand  thermal  effects  can  be  neglected  when  the 
energy  spread  in  the  beam  S,  E  «  Eq/  N,  where  E0  is  the  bulk  beam  energy  and  /V  is  the  number  of 
wiggler  periods  in  the  interaction  region.  In  this  regime  the  entire  electron  beam  participates  in  the 
interaction;  however,  collective  (i.e.,  electrostatic)  effects  are  important  only  when  the  fluctuating 
space-charge  potential  is  comparable  to  the  ponderomotive  potential.  The  interaction  in  this  collective 
regime  is  referred  to  as  stimulated  Raman  scattering,  and  describes  the  coupling  of  a  negative  energy 
space-charge  wave  and  a  positive  energy  electromagnetic  wave  through  the  presence  of  the  wiggler.  On 
the  other  hand,  when  thermal  effects  are  important,  the  radiation  is  resonant  with  only  a  small  fraction 
of  the  beam  and  the  process  is  termed  stimulated  Compton  scattering.  We  shall  be  concerned  in  this 
paper  with  the  cold  beam  limit,  and  deal  with  both  the  single-particle  and  collective  regimes. 
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An  additional  factor  in  the  interaction  is  introduced  by  practical  limitations  in  the  propagation  of 
intense  electron  beams;  specifically,  that  an  axial  guide  field  is  required  to  collimate  intense  beams  in 
the  transverse  direction.  This  is  of  primary  importance  to  free  electron  laser  experiments  operating  at 
millimeter  wavelengths1-5  which  employ  relatively  high  current  (>  1  kA)  and  low  energy  (~1  MeV) 
electron  beams.  In  contrast,  an  axial  guide  field  is  not  a  practical  necessity  for  infrared  free  electron 
laser  experiments6-8  which  typically  operate  at  much  lower  ambient  currents  (~  1  A)  but  higher  ener¬ 
gies  (- —  50  MeV).  However,  the  effect  of  the  axial  guide  field  on  the  free  electron  laser  mechanism 
may  be  relevant  in  the  latter  case  as  well  because  of  enhancements  in  the  gain  which  may  result  when  a 

guide  field  is  present. 

The  effects  of  an  axial  guide  field  have  been  treated  from  the  standpoints  of  both  a  fluid9-13  and  a 
kinetic14-16  theory.  It  is  our  purpose  in  this  work  to  treat  the  question  of  the  free  electron  laser  insta¬ 
bility  in  the  presence  of  an  axial  guide  field  in  both  the  tenuous  (i.e.,  single-particle  regime)  and  dense 
(i.e.,  Raman  regime)  beam  limits  by  means  of  a  solution  of  the  Vlasov-Maxwell  equations.  In  the 
interest  of  analytic  tractability  a  cold-beam  approximation  is  imposed,  and  the  contributions  due  to 
cyclotron  mode  interactions  are  included.  The  organization  of  the  paper  is  as  follows.  In  Sec.  II.  we 
develop  the  fluctuating  source  currents  by  solution  of  the  Vlasov  equation.  The  unperturbed  orbits  are 
assumed  to  be  constant  axial  velocity  (helical)  trajectories,17  18  and  the  source  currents  are  found  for  a 
general  equilibrium  distribution.  The  general  dispersion  equation  is  obtained  in  Sec.  Ill,  and  solved  in 
several  analytically  accessible  regimes  for  a  cold  beam  limit.  A  detailed  numerical  solution  is  presented 
in  Sec.  IV  for  a  wide  range  of  operating  parameters.  A  summary  and  discussion  appears  in  Sec.  V. 

II.  THE  SOURCE  CURRENT 

In  this  section,  we  derive  the  fluctuating  source  current  by  means  of  solution  of  the  linearized 
Vlasov  equation.  The  physical  configuration  we  consider  is  that  of  a  relativistic  electron  beam  propagat¬ 
ing  through  an  ambient  magnetic  field  composed  of  a  periodic  helical  wiggler  field  and  a  uniform  axial 
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B  -  ZJ0e.  +  Bw(/-.  z).  (1) 

where  the  wiggler  field  is  assumed  to  be  generated  with  a  bifilar  helix  and  is  derived  from  a  vector 

potential  of  the  form19 

A,- - fMtV  cos  -  l\(kwr)  sin  (0- Jtwz)e#j  (2) 

in  cylindrical  coordinates.  In  Eq.  (2),  Bw  and  kw  (=  2ir/Aw,  where  is  the  wiggler  period)  are 

assumed  to  be  constant  and  I\  and  l[  are  the  modified  Bessel  function  of  the  first  kind  and  its  deriva¬ 

tive  respectively.  Since  for  most  free  electron  laser  experiments  the  initial  beam  radius  is  a  small  frac¬ 
tion  of  the  wiggler  period,  we  shall  expand  in  powers  of  £wrand  write 

B„  =*  Bw(hx  cos  kwz- t-e„  sin  kwz).  (3) 

The  single-particle  orbits  in  these  combined  fields  have  been  amply  discussed  in  the  literature,17  18 
and  will  not  be  discussed  in-depth  here.  We  shall  restrict  consideration  to  orbits  which  are  given 
approximately  by  stable  helical  trajectories,  and  write15- 18 

px  =-  ymvw  cos  kwz  +  Px  cos  ft  o'  -  P,  sin  ft0r. 

p,»ymv,  sin  kwz  +  Px  sin  Q0t  +  Pv  cos  ft0r.  (4) 

p.  -  ymvn  -  P„[PX  cos  ( kKz  -  ft0f)  -  P>  sin  (kwz  -  ft0r)L 
where  vw  =  ft  W\J  (ft0  -  Awv„)  and  v„  are  the  transverse  and  axial  velocities  corresponding  to  the  heli¬ 
cal  trajectories,  ft0.„  =  \eBnJymc  |,  y  =  (1  -  v2/c2)-1/2,  Bw  =  vw/vn  is  the  pump  strength  parameter, 
and  Px  and  Py  are  constants  which  correspond  to  the  canonical  momenta  in  the  limit  as  B0  —  0.  Equa¬ 
tions  (4)  are  valid  as  long  as  1/3 WPX,V\  «  ly/nvul  and  require  that  \w  and  vM  are  related  via 

W  +  v2  -  (l  -  y-2)c2,  (5) 

which  constitutes  a  quartic  equation  for  vt, .  Equation  (5)  describes  at  most  four  distinct  classes  of  tra¬ 
jectories  of  which  one  is  characterized  by  motion  antiparallel  to  B0  and  will  be  ignored.  Of  the  remain¬ 
ing  trajectories  having  motion  parallel  to  B0,  we  restrict  consideration  to  those  which  are  stable,17 18  i.e., 
for  which 

(A\,v„  — ftoM/fc^v,,-  (l+j32)ft0l  >  0.  (6) 


3 


FREUND,  SPRANGLE,  DILLENBURG,  da  JORNADA,  SCHNEIDER.  AND  LIBERMAN 


The  typical  dependence  of  the  axial  velocity  on  B0  is  shown  in  Fig.  1,  in  which  we  plot 
j8||(s  v,,/c)  versus /30  (=  fl^kwc)  for  a  kinetic  energy  of  1.5  MeV  and  wiggler  amplitude  and  period 
such  that  fi  J kwc  —  .05.  There  are  two  classes  of  stable  orbits.  One  class  (referred  to  as  group  I)  of 
orbits  is  characterized  by  high  axial  velocities  and  has  kwvu  >  (1  +  /3»)fi0.  For  these  trajectories  the 
axial  velocity  decreases  monotonically  with  increasing  B0  until  vM  —  (1  +  0  1)0  qC  which  is  the  orbital 
stability  boundary.  For  the  parameters  shown,  this  occurs  for  0O  =  .765.  The  second  class  of  orbits 
(group  II)  is  characterized  by  an  axial  velocity  which  increases  monotonically  from  zero  with  increasing 
B0%  and  high  axial  velocities  are,  typically,  found  only  when  /30  >  1.  In  contrast  to  group  I  trajectories, 
the  stability  criterion  (6)  is  satisfied  for  these  orbits  because  kw\, ,  <  ft0.  It  will  be  shown  in  Sec.  Ill 
that  the  relative  magnitude  of  kwvu  and  fi0  for  group  I  and  II  orbits  has  important  consequences  on  the 
excited  spectrum.  Finally,  it  is  important  to  note  that  £»v.,  fl0  f°r  either  class  of  orbits.  Were  the 
equality  to  hold  it  would  imply  an  infinite  vw  which  would  violate  the  conservation  of  energy 
represented  by  Eq.  (5). 


The  source  current  is  obtained  from  the  velocity  moments  of  the  perturbed  distribution  function 
8/6(z.  p.  r)  —  fb(z.  p.  t )  -  Fb(Px,  Pv,  p),  where  f„  is  the  complete  distribution,  Fb  is  the  equilibrium 
distribution,  and  S  fb  is  assumed  to  be  first  order  in  the  radiation  fields.  The  formal  solution  of  the 
Vlasov  equation  to  this  order  is 


S/4(z,  p,  r(z))  -  e 


r  dz' 

8E(z\  r (z ))  +  -  v(z’)  x  8B(zr.  Mz')) 

v.(z') 

c 

0Fb 

dp(z')  ’ 


(7) 


where  the  solution  is  parametrized  in  terms  of  the  axial  position  relative  to  the  start  of  the  interaction 
region  (at  z  ”  0)  and  Mr)  *  t0  +  dz'/wAz')  is  the  sum  of  the  time  required  for  an  electron  to 
transverse  the  distance  z  and  the  entry  time  t0. 


We  assume  plane  wave  solutions  of  the  form  exp(-/u>f)  and  choose  to  work  with  the  scalar  and 
vector  potentials 


8<6  (z.  t ) 


—  8<b(z)  exp(—  io»f)  +  c.c.. 
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8 A(z.  t)  —  -j  8  A(z)  exp(- i(ot)  +  C.C.. 

where  e,  •  8  A(z)  **  0.  In  addition,  since  the  assumption  of  small  Px  and  Py  are  implicit  in  the  analysis, 
we  adopt  an  equilibrium  distribution  of  the  form 

Fb(Px.  Py.  p)  -  nb8(Px)8(Py)Gb(p).  (8) 

where  nb  is  the  number  density  of  the  beam,  and  Gb(p)  is  an  arbitrary  function  subject  to  the  normali¬ 
zation  condition  J*Q  dppGp(p)/p:  *  1.  Observe  that  the  choice  of  a  distribution  of  this  form  (8) 
confines  the  equilibrium  trajectories  to  be  those  constant  axial  velocity  (helical)  orbits  described  by  (5). 
The  interested  reader  is  referred  to  Freund  et  ai.!5  for  a  detailed  derivation  of  the  perturbed  distribution 

and  source  current.  We  confine  ourselves  here  to  the  final  result.  With  respect  to  the  basis  ex  =*  ~ 
(ev  ±  ie„),  the  source  current  can  be  written  in  the  form 

8  J(z,  r)  —  (8y+(z)e+  +  8/_(z)e_  +  8/.(z)erl  exp (—/«/)  +  c.c..  (9) 


where 


SJ±  (z)  -  fQ  dp  -y-  exp(T/fl0/(z))  2  +  — y~  £>*  +  exp(±t'ft0/(z))  —  D- 


a 

+  i  D+  + 

a 

.  a 

D--  D.~- 

Gb{p ) 

dPx 

BP, 

BPX 

'  a  py 

■  bp 

.  (ll) 


In  Eqs.  (10)  and  (11).  <uft  is  the  beam  plasma  frequency,  =  px^  ipy , 

D±  =  -  exp[±/ft0' (->H8 -‘i±(‘)  ~  8 .4  *  (0)  exp [/(« T  fl0)r (z.  0)] 

^/n0  f  .4 ±(z')  H0)t (z.  z')]\, 

J0  v.(z) 


D,  =  J*  dz  exp  [iuir  (z,  z')]  -3,  8 i(z')  +  —  ~8J  +  (z')  -f  (13) 

P  J°  ‘  c  p.  p.  j| 
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where  r(z.  z')  =  r(z)  -  r(z’).  Finally,  we  assume  the  spatial  dependence  of  the  fields  to  be 


8  A  ±  (z)  —  8,4  s  (0)  exp  Uk±z) 


8 <Mz)  -  8^(0)  exp  (ikz). 


As  a  result,  the  source  currents  can  be  expressed  as 


8/±(z)~-  ^S^fzlL/A-  -7.-77^  aiS±ik±.m) 
me  \  pz  <u  +  tip  -  XiV,,  // 


-  a*  -r—  8<i  (z)<?  ' '\a{k.  <o), 


87.  (z)  -  -  ^8<£(z)*(A,  tu)  -  -3— ^  aH.fs/i^(z)  exp  UkHlz)cr+(k+.  to) 
8ir  on-  l 

+  8/t_(z)  exp  (—ikw:)<rAk-,  <u)|, 

J,o° 

Gb(p)(— ), 


x(fc.  to)  —  ~  f  dp  - - 

A  •'0  cu  — 

,,  ,  _  “t  f"  ,  *WV'I 

*.<*.»>  =  —  J„  4-  t-vl;  -  ri„ 

5±(^±,aj)=/  dp—P—  Gb(p)  — - 2— j[l 

n  y />n  Uvv,,-  ftp)  1 


,  _  <“l  r°°  j  /w  dGf> 

“  ~  zt  Jo  P  w  -  Av„  dp 

>l_  C°°  ,  *»rV'l _ AW 

A  ^  A»v i|  —  ftp  y(cij  —  AvH) 

,  .  c2A2  r. 


Op  (P  ^  fl q  ±  A„.  v,| _ 01  T  Ar^v„ 

A^v,,  -  12 0  aiTUp  -  Av„  to  kw Vq  -  Asv,, 

,  r°°  .  A2v,,  m  d 


■ST* 


0  y2(AwvM  -  fi0)2  u,*kwv;t-  k±v„  dp 


oj  i  r 00 

rr±(k±,  to)  =  —  I  dp 
(j)  ^  0 


y(/rH.v„  -  n0)  to*kMvl{  -  A±v, | 


yVn  (A^v,,  -  n0)2 1  <0 T Arwv.,  -  k±\,,  tu^Op  -  Atv., 


and  it  has  been  assumed  that  Gb(p )  -=  0  for  p  —  0,  00.  Observe  that  Eqs.  ( 14)-<  19)  are  equivalent  to 
the  result  found  by  Sprangle  and  Smith20  in  the  limit  as  B0  —  0. 
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III.  THE  DISPERSION  EQUATION 

The  dispersion  equation  is  obtained  by  substitution  of  the  source  currents  [Eqs.  (14)- (19)1  into 
Maxwell’s  equations 


a.2  + 


&y|±(z)-  -  ~bJ±(z). 


8, S<Mz)  -  —  8/.(z). 

0) 

This  results  in  a  set  of  three  coupled  equations  for  the  initial  field  amplitudes 


(20) 


A  1  • 

A-Uc  kw,  <o)  8A±(0)  +  -r  aw  — r-  *a(/c,  <o )  (0)  =■  0, 

2  o> 

k  A 

ll  +  x  (k.  <*»)!  8<£  (0)  +  [8^  +  (0)<t  +  (^  —  k„,  a>) 

+  si_(0 )(T_(k  +  kw,  co)]  -  0. 

where  the  wavelength  matching  conditions  imply  that  k±  -  k^kw,  and 


(21) 


\  +  (k±.  w)  s  l  - 


c2kl  -2 


_  'IL/JL. 


(O  -  k±  Vi  I 


uj-  a )-  \y Pz  uTfl0-(tv I,  /  2a>-’ 


-  t  ~  al  S±(k ±.  to) 


(22) 


,2  ...2 

describes  the  dispersion  properties  of  the  pure  electromagnetic  modes  in  the  combined  wiggler  and  axial 
guide  fields.  The  dispersion  equation  itself  is  found  by  setting  the  determinant  of  the  matrix  of 
coefficients  of  Eqs.  (21)  to  zero,  which  yields 


,  ,  .  «l  kw  c2k:  ,,  . 

1  +  \(k  co)  -  — - - - —  \a(k.  <o) 

2  kat- 


cr+{k  —  kK.  co)  <r_(k  +  /fB.  co) 


(23) 


\+(k  —  kv.  co)  A_(Ar  +  kw.  co) 

It  is  evident  that  Eq.  (23)  describes  the  coupling  between  the  electrostatic  beam  mode  with  each  of  the 


electromagnetic  modes. 


In  the  interest  of  analytic  tractability,  we  shall  now  assume  that  the  electron  beam  is  sufficiently 
cold  that  a  monoenergetic  distribution  of  the  form 

<j*(p)  “  —  S(p  -  pn)  (24) 

P 

can  be  employed.  As  mentioned  previously,  this  is  generally  valid  as  long  as  the  momentum  spread 
\p  «  PrJ  N.  Combination  of  (23)  and  (24)  yields  a  dispersion  equation  of  the  form 


7 


FREUND.  SPRANGLE.  DILLENBURG.  da  JORNADA.  SCHNEIDER.  AND  LIBERMAN 


,  ,  ,  Pi  u>i  a*(A  —  kw.  w)  «_(A  +  ku,m) 

(«  -  kv.y-KW-  -  —  — n — r — r +  "-71-  r~~  r  (2S) 

2  y  e  +  \k  —  A*,  co)  «_(A  +  A„  ,<o ) 

to  second  order  in  the  wiggier  amplitude,  where  all  orbit  quantities  (i.e.,  v., ,  y..  fi„.  etc.)  are  computed 
using  Po. 


11-  * 

K  V,f  =  - -  <P. 

yy- 


« ~(k±.  co)  =  1  — 


co j(co  ~  k  a  V11) 

y(co2  —  kic2)  ((uTfli)-  Asv,, ) 


<x±(k  T  A...  co)  =  — - — — ry  £i,2(co2  -  A2c2  -  u>l/y) 

oi  —  (A:  +  Aj-c2 


+  —  /3,2(co  -  A' v, )  2  ((1  -  'I'2)  (At  -  coV|,/c)2  +  y.  2(1  -  4>)  (co2  -  A‘r2)l 

y 

Co  l  fl  n  CO  f  A*.  Vi.  (jj 

+  — - 2— -  /3„<At  -  cov„/c)  'P  - — - — - ■■  "  — 

you  A^.v,,  —  fl0  co  —  Av„  uj^CIq—  (A— A„.  )vM 


3co  T  AK.v., 


co  T  n0  i  AH  v,. 


**0  r,  .  ,2  1  21  "  **0 

+  T - 77-  [(«>  -  Ac v;, > 2  -  *2v„2]  - - - - - - 

A^v,  —  ft0  co  —  A vu  <o  T  Q  0  —  (A  T  A„)v 


—  y.  2(1  —  <P)  [(co  —  Av,|)2  —  «2Vi2] 


In  addition,  we  have  defined 


co  <o  -  2A Vu  2ft0 

co  —  Av„  co  —  Ay.,  A*  v,  —  n, 


<t>  =  1 


(1  +  p 2)  fi0  -  Awvn  ’ 


n0  (1  +  /32.)  Aw Vu  -  co 

vj,  =  1 - - - — - — - .  (30) 

fiw  ( Ac  -  CO  v,,/ c )  ( l  +  j3 2 )  Q  0  -  A„. V, 

In  Eq.  (25),  c.lAT  A*,  co)  describes  the  circularly  polarized  electromagnetic  modes  in  the  absence  of 


the  wiggier,  and  the  left-hand-side  describes  the  electrostatic  beam  modes.  Observe,  however,  that  the 


presence  of  the  wiggier  modifies  the  natural  electrostatic  response  frequency  by  a  factor  of  <P  defined  in 


The  wiggier  field,  therefore,  provides  a  coupling  between  the  space-charge  wave  and  either 

polarization  state  of  the  electromagnetic  wave.  We  choose,  without  loss  of  generality,  to  focus  on  the 

coupling  with  the  3/1*  mode.  As  a  consequence,  we  shall  assume  that  l«,(A-  A„.  co)|  «  !e_(A  + 
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A'h,  <*»)|  and  neglect  the  term  in  e_  l(k  4-  kw,  u)  in  (25).  If  we  assume  in  addition  that  u>l/yoj2  «  1, 
then  the  dispersion  equation  can  ho  cast  into  the  substantially  simpler  form 


[(<*  -  /cvi()2  -  k 2Vi2l  to2  -  klc2  - 


-  A;+V||) 

y  (<u  —  ft  o  “  *+Vn) 


& b 


(uhtu  -  <C+V||) 


ZL  a  2  ^2  _  k2c2  _  _ 

—  —T~  Pit  to  —  K  C  —  -  ~  - 

2  y  y(tu-n0  -  k+v, ,) 


4-  y.~2ti)2(l  —  <l>)  +  o*n 


to  —  Arvii 

0  to  -  n0  -  k  +  Vi | 


where  we  have  written  k+  -  k  -  it*  for  simplicity.  Peak  gain  in  (31)  can  be  expected  to  occur  near  the 
intersections  of  the  electrostatic  and  electromagnetic  dispersion  curves.  A  schematic  representation  of 
the  dispersion  relation  is  shown  in  Fig.  2  for  fl0  >  to4/y1/2.  Evidently,  high  frequency  interactions 
with  the  electrostatic  beam  mode  can  occur  only  in  the  positive  co  and  k+  quadrant  when  v,  >  0,  and 
we  shall  restrict  the  analysis  to  this  regime.  It  should  be  observed  that  the  slope  of  the  space-charge 
modes  (to  =*  k+vn  +  fc„Vi,  ±  kv„)  is  identical  to  the  pure  cyclotron  mode  (a>  —  A:^ v,  +  tl0),  and  the 
relative  magnitudes  of  the  k+  -  0  intercepts  of  these  curves  determine  which  branch  of  the  electromag¬ 
netic  dispersion  curve  participates  in  the  interaction.  This  point  will  be  discussed  in  more  depth  at  a 
later  stage  of  the  analysis. 

The  dispersion  equation  represented  by  (31)  is  a  fifth  degree  polynomial  in  k.  If  we  make  the  res¬ 
triction  that  k+  >  0,  then  (31)  can  be  reduced  to  the  following  quartic  equation 

8  2 

(k  —  — —  k)  (k  —  —  +  k)  (k  —  kw—  KJ  {k  -  kw  -  KJ  =  -  ~~~  S2k?  ~r  8~' 

V|,  Vt,  4  Kc 


x  -^-4>k-kw-  ^  -Sll  k  _  «.  , 

y-V.,  V.,  V;,  Vn 


where  £  =  u>blyU2ckn  is  the  beam  strength  parameter,  c2K 2  =  a*2  -  o>A2/y, 


_  1  ^  ,  w  —  Q  o  I  _  1 


K±  =  y  K  + 


±  x  \K2  +  2t2k 


2  4.  1(2^2  ill 


and  IK  S  K  -  (<o  -  fi0)/vi.  It  is  clear  that  the  nature  of  the  interaction  is  strongly  dependent  upon 
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the  sign  of  <t>,  which  affects  the  natural  electrostatic  response  frequency  of  the  plasma  as  well  as  mediat¬ 
ing  the  ponderomotive  force.2  In  the  limit  as  B0  — •  0,  <t>  approaches  unity  and  (32)  reduces  to  the 
well-known  results  in  the  limit  of  zero  axial  field.20-22  The  behavior  of  d>  when  the  guide  field  is  finite, 
however,  is  strongly  dependent  on  the  type  of  orbit  under  consideration.  The  variation  of  <t>  with/30  is 
shown  in  Fig.  3  for  the  parameters  used  to  generate  the  orbits  in  Fig.  1.  As  shown  in  the  figure,  <1>  ^  1 
for  group  I  orbits  and  contains  a  singularity  when  fl 0(  1  +j82)—  /c„.vM  which  is  the  orbital  stability 
boundary.  In  the  case  of  group  II  orbits  the  magnitude  of  <l>  is,  typically,  less  than  or  of  the  order  of 
unity,  but  <!>  is  negative  for  axial  guide  fields  less  than  a  critical  magnitude  given  by 
[1  +/32(l-y2)]fi0“  k  wv,|.  As  B0  increases  beyond  this  critical  value  <f>  approaches  unity;  however, 
it  should  be  noted  that  *.  decreases  monotonically  with  increasing  B0  in  this  regime.  In  either  case  in 
which  <l>  >  0,  the  interaction  is  basically  one  in  which  a  positive  energy  electromagnetic  wave  is  cou¬ 
pled  to  a  negative  energy  space-charge  wave  by  the  action  of  the  wiggler.  Neither  wave  is  intrinsically 
unstable  and  growth  occurs  when  the  wiggler  amplitude  is  above  threshold.  However,  when  <1>  is  less 
than  zero  k  is  imaginary  and  the  space-charge  waves,  themselves,  comprise  a  complex  conjugate  pair 
u>  =»  Jcv„  ±  /UivM  one  of  which  is  unstable.  As  a  result,  we  shall  distinguish  between  these  two  possi¬ 
bilities  and  treat  the  solution  to  the  dispersion  equation  when  <l>  is  positive  and  negative  separately. 


A.  d>  >  0 


In  this  regime  we  observe  that  the  orbital  stability  criterion  implies  that 

12/31 


k*v„  >  n0 


l  + 


5. 


A 


for  group  I  orbits,  and  the  requirement  that  <t>  be  positive  leads  to  the  condition  that 


<  n0 


1  -  (y;-  1),/3 


A 

2/3 

Bo 

(34) 


(35) 


for  group  II  orbits.  As  a  consequence,  the  intersection  between  the  space-charge  and  electromagnetic 
modes  occurs  at  frequencies  greater  (less)  than  ft0(l  —  >3,, )~ 1  for  group  I  (II)  orbits  when  *v„  « 
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kwv.|  -  n0|.  This  condition  is  satisfied  as  long  as  <t>  <  y.-f  '(BJ B0)13  and,  since  $  «  1  is  impli¬ 
citly  assumed  in  order  to  neglect  self-field  effects,  is  a  relatively  weak  constraint.  A  schematic  represen¬ 
tation  of  the  interaction  is  shown  in  Figs.  4a  and  4b  for  group  I  and  II  orbits  in  which  we  plot  <u  versus 
and  the  dotted  line  represents  complex  conjugate  roots.  It  is  the  lower  (i.e.,  u>  “  Xc  vM  -kv„) 
space-charge  mode  which  produces  the  active  coupling  and  wave  growth  since  this  is  the  negative 
energy  mode.  We  note  that  from  (34)  and  (35),  the  intersections  are  not  close  to  the  cyclotron  line 
and  occur  approximately  for  &>  =*  ck+  =  (k+  +  kw)\».  This  is  the  well-known  free  electron  laser  reso¬ 
nance  at  k+  =  2y.2&j3n. 

If  the  beam  strength  parameter  is  sufficiently  small  that  f  «  y.  (BJB0)Vi  and  yABj  B0)V} 
<l> l/2,  then  the  cyclotron  resonance  effects  can  be  neglected  and  the  dispersion  equation  reduces  to  the 
cubic 

a  2 

8k{8k  +  2k)(8 k  -  \k)  =  -  ^2/cj0,r '<!>-— -  (36) 

4  y;v.i 

where  we  have  chosen  8k  =  k  -  co/v„  -  k,  and  Ik  =  kw  +  K  -  w/v,,  -  k  is  the  frequency  mismatch 
parameter.  Equation  (36)  reduces  to  the  result  found  by  Sprangle  and  Smith20  in  the  limit  as  B0  — -  0, 
and  corresponds  to  the  limit  discussed  by  Friedland  and  Fruchtman.13 


The  single-particle  or  "strong-pump"  regime  occurs  when  |5  A I  »  |2tc|.  In  this  regime,  (36)  can 
be  approximated  by 


8 k2(8k  -  bk)  =  -  ^rt2kJ3 ,r‘<t> 

4  y/v„ 

Peak  growth  occurs  when  Ik  =  0,  for  which  the  complex  roots  are 


(37) 


<Sk)ma*  =  y  (v3  ±  f) 


Mi 


kw. 


(38) 


Therefore,  self-consistency  imposes  the  requirement  that 


x  «  -jy/3 lylP'i  kw 


(39) 


in  order  for  Eq.  (37)  to  be  valid.  It  is  important  to  recognize  that  in  this  regime  the  coupling  between 
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the  electrostatic  beam  mode  and  the  6 .4  ^  mode  is  relatively  unimportant,  and  the  strength  of  the  pump 
and  ponderomotive  potential  completely  dominate  the  interaction.  Because  of  this,  (37)  can  be 
recovered  in  a  much  more  direct  manner  by  ignoring  the  space-charge  potential  in  (21)  and  setting 

A Ak  -  kw.o>)  -  0. 


B.  <I>  <  0 


As  mentioned  previously,  the  space-charge  waves  are  intrinsically  unstable  in  this  regime  and  a 
relatively  broadbanded  spectrum  of  excited  electrostatic/electromagnetic  waves  is  expected  to  occur.  In 
addition;  since  we  are  dealing  with  group  II  orbits 

B  2/3 

fto  1  ~  (y]  ~  1)1/3  -r—  <  **.v„  <  n0.  (44) 

O0 

and  the  possibility  of  a  cyclotron  mode  interaction  exists  when  y?  -  1  «  t  BJ Bw)2  . 

In  the  high  frequency  <a  »  (l-j3)_1n0  and  high  k  »  k„  regime,  the  space-charge  waves  are 
characterized  by  frequencies  at  «  ck *,  and  the  cyclotron  resonance  contributions  can  be  ignored  when 
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£  «  |<t> |l/:/3„y_.  '(-y2-  1)1/3 

ck„ 


Bw 


B  o 


2/3 


(45) 


In  this  regime  we  recover  the  cubic  dispersion  Eq.  (36)  found  earlier  (when  4>  >  0).  For  frequencies 
such  that  \lk\  «  8k  the  free-electron  laser  solutions  are  obtained  [see,  for  example,  Eqs.  (37)  and 
(40)  for  the  strong-pump  and  space-charge  dominated  regimes);  however,  for  high  frequency  space- 
charge  modes  I A  A:  |  =  cu(l  —  /3„ )/  v^  »  \8k\.  For  these  waves,  the  dispersion  equation  can  be  approx¬ 
imated  by  the  quadratic 


2 


k2  -  2  —k  +  +  k  I 


Pi 

1  -  +pn) 


—  0. 


vn  v,r 

which  has  been  expressed  in  terms  of  k  rather  than  8  k  for  convenience.  The  solution 


(46) 


*  —  —  ±  / 1*  I  -v  / 1  -  ^-PnyHl  +0„).  (47) 

v,j  V  4  - 

is  obviously  a  modified  space-charge  wave.  It  should  be  observed  that  the  presence  of  the  wiggler  acts 
as  a  stabilizing  influence,  and  for  sufficiently  strong  pumps  (i.e.,  P lP»yHl  +  P»)  >  4)  the  mode  is 
stable. 


In  the  opposite  limit  in  which  <u  «  ( 1  —  >3U )~ 1 Q  0,  it  is  more  difficult  to  satisfy  condition  (45), 
and  the  cyclotron  resonance  is  of  greater  importance.  This  regime  will  be  discussed  in  depth  in  Section 
IV  in  the  context  of  a  complete  numerical  solution  of  the  dispersion  equation. 


IV.  NUMERICAL  ANALYSIS 

In  this  section,  we  conduct  an  extensive  numerical  analysis  of  the  quartic  dispersion  Eq.  (32)  for  a 
beam  characterized  by  -y  *■*  3.94  (1.5  MeV)  and  {  —  .1,  and  a  wiggler  such  that  Q  J kKc  —  .05.  Our 
procedure  is  to  solve  (32)  for  a  wide  range  of  axial  field  strengths  in  the  vicinity  of  fl0  =  kwc  by  self- 
consistently  calculating  v.,  for  each  case  and  for  both  types  of  stable  trajectory.  As  in  Section  III,  we 
distinguish  between  the  regimes  for  which  <1>  is  positive  and  negative. 
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In  the  limits  in  which  4>  is  positive,  the  dispersion  properties  of  (32)  are  qualitatively  represented 
in  Fig.  4  for  orbits  in  groups  I  and  II.  The  growth  rates  (Im£//c„)  are  plotted  versus  frequency  for 
several  appropriate  values  o(  0O  in  Figs.  5  and  6  for  the  two  types  of  trajectory.  Instability  is  found  to 
occur  for  o>  =  (1  -  0u)~lkwV;  which  corresponds  to  the  well-known  free-electron  laser  resonance  con¬ 
dition.  The  observation  of  decreasing  (increasing)  frequencies  of  instability  with  increasing  axial  field 
for  group  I  (II)  orbits  can,  therefore,  be  explained  by  examination  of  Fig.  1  in  which  it  is  seen  that  v, 
decreases  (increases)  with  increasing  0O- 

The  behavior  of  the  peak  growth  rates  for  the  two  classes  of  orbits  can  also  be  readily  explained. 
As  shown  in  Fig.  7,  the  peak  growth  rate  for  group  I  orbits  is  a  monotonically  increasing  function  of  Bn 
up  to  the  singularity  at  the  orbital  stability  boundary  (at/30  —  765  for  this  choice  of  parameters).  This 
behavior  is  due  to  increases  in  both  0  w  and  d>  with  the  axial  field  which  results  in  increases  in  the 
effects  of  both  the  electrostatic/electromagnetic  coupling  and  the  ponderomotive  potential.  Observe, 
however,  that  the  singularity  is  due  solely  to  the  character  of  d>  since  0  „  is  everywhere  finite.  The  scal¬ 
ing  of  the  maximum  growth  rate  with0o  is  also  shown  in  Fig.  7  for  group  II  trajectories.  In  this  case, 
however,  <t>  is  bounded  by  unity  and  increases  monotonically  from  zero  (at0o  =  1.25)  with  increasing 
axial  fields.  In  addition,  0W  decreases  monotonically  to  zero  with  increasing  fl0  for  group  II  orbits  since 

lim  0W“  BJ  Bn.  As  a  consequence,  the  peak  growth  can  be  expected  to  initially  increase  from  zero  at 
Bn~' » 

0o  =  1.25,  and  to  decrease  again  slowly  to  zero  as  the  axial  field  becomes  large.  This  behavior  for  the 
growth  rates  of  each  class  of  orbit  is  in  qualitative  agreement  with  that  found  previously  in  the  context 
of  a  low  gain  theory.15  Finally,  since  it  is  our  intention  to  treat  the  collective  regime,  the  parameters 
have  been  chosen  to  correspond  to  the  space-charge  dominated  limit  discussed  in  Section  III  and  the 
numerical  results  for  the  peak  gain  can  be  recovered  from  Eq.  (40)  to  within  an  error  of  a  few  percent. 

Some  care  must  be  taken  in  the  characterization  of  the  dispersion  properties  of  (32)  when  <  0. 
For  low  axial  fields  (j3 0  <  .7  for  the  parameters  under  consideration)  both  v„  and  14*1  are  low  and  con¬ 
dition  (43)  is  not  well-satisfied  when  <u  <  (1  -  0,)-lno.  This  regime  is  illustrated  in  Figs.  8  and  9  in 
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which  we  plot  oil  ckw  versus  R ekjkw  and  lmk/kw  versus  a >/ckw  respectively  for0o  "  -5.  The  dashed 
line  in  Fig.  8  corresponds  to  complex  conjugate  roots.  Evidently,  two  instability  regimes  exist.  At  high 
frequencies,  the  modified  space-charge  wave  discussed  in  Section  III  B  is  obtained  with  an  asymptotic 
value  (i.e..  high-*)  limit  of  lm  klkw)  which  agrees  to  within  1%  of  the  predicted  value  in  Eq.  145). 
The  instability  found  at  lower  frequencies  is  difficult  to  treat  analytically,  and  corresponds  to  a  modified 
cyclotron  mode. 

For  higher  axial  fields  1/3 0  >  .7),  the  character  of  the  unstable  modes  is  altered.  As  shown  in  Fig. 
10,  in  which  we  plot  oil ckw  versus  RskJkK  for/30  -  1,  there  are  still  two  unstable  regimes.  While  the 
higher  frequency  regimes  corresponds  to  the  modified  space-charge  mode  in  this  case  as  well,  the  lower 
frequency  instability  requires  some  discussion.  The  growth  rate  in  this  regime  is  plotted  as  a  function 
of  frequency  in  Fig.  11  for£0  =■  .8,  1,  and  1.2.  In  each  case,  the  peak  in  lmk/kw  observed  for  the 
lower  frequency  instability  corresponds  to  the  region  shown  in  Fig.  10  in  which  the  real  part  of  the 
complex  modes  (dashed  line)  exceeds  that  of  the  waves  which  are  purely  real.  Near  the  peak  Ik  =  0, 
and  the  instability  which  results  is  a  largely  electromagnetic  free  electron  laser  interaction  which  (for 
our  choice  of  parameters)  agrees  to  within  a  few  percent  of  the  analytic  expression  for  the  peak  growth 
rate  (40)  in  the  collective,  or  stimulated  Raman  scattering,  limit.  As  the  frequency  decreases.  I A  A I 
increases  and  the  character  of  the  instability  becomes  increasingly  electrostatic  and  we  find  an  unstable 
modified  space-charge  wave  for  frequencies  o>  >  ft0. 

V.  SUMMARY  AND  DISCUSSION 

In  this  paper  we  have  analyzed  the  linear  growth  rate  in  both  the  single-particle  and  collective 
regimes  of  operation  of  a  free  electron  laser  configuration  which  contains  a  uniform  axial  guide  field. 
The  technique  employed  consists,  essentially,  of  a  Vlasov  theory  of  the  perturbations  about  constant-v,. 
helical  trajectories,  and  includes  the  effects  of  both  stimulated  Raman  scattering  with  electrostatic  beam 
modes  and  the  effect  of  the  ponderomotive  potential  due  to  the  excited  radiation.  Analytic  expressions 
for  the  growth  rate  in  these  two  regimes  are  found  and  a  comparison  is  made  with  a  numerical  solution 
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of  the  full  dispersion  equation.  Substantial  agreement  is  found  between  the  analytical  and  numerical 
treatments.  In  addition,  substantial  qualitative  agreement  also  exists  between  the  present  work,  which 
is  in  the  high  gain  regime,  and  the  small  signal  theories.1  '• 15 

The  results  indicate  that  substantial  enhancements  in  the  growth  rate  of  the  free  electron  laser 
instability  may  be  obtained  by  the  inclusion  of  an  axial  guide  field  in  which  ft0  —  as  a  result  of 
increases  in  both  the  transverse  velocity  and  the  ponderomotive  potential.  It  should  be  pointed  out, 
however,  that  the  presence  of  the  guide  field  also  gives  rise  to  instability  in  the  electrostatic  beam  mode 
for  group  II  orbits  when  d>  <  0.  Since  this  may  have  a  degrading  effect  on  beam  quality  and,  in  turn, 
on  the  operating  efficiency,  further  study  of  this  question  is  required. 

A  cyclotron  mode  interaction  in  the  collective  regime  has  been  found  only  where  the  axial  velo¬ 
city  and  !<t>  I  are  relatively  small.  However,  two  factors  should  be  noted.  The  first  is  that  our  choice  of 
equilibrium  orbits  (i.e.,  with  Px  -  Pv  -  0)  has  intrinsically  ignored  the  random  (or  Larmor)  component 
of  the  transverse  velocity  which  can  be  expected  to  be  the  primary  source  of  a  gyrotron  type  of  instabil¬ 
ity.  The  second  is  our  choice  of  an  idealized  magnetic  field  structure  in  which  transverse  gradients  in 
the  wiggler  field  have  been  neglected,  which  is  valid  only  for  orbits  in  which  k„r  «  1.  Since  small  v, 
implies  large  transverse  velocities  and  excursions  from  the  axis  of  symmetry,  the  choice  of  an  ideal 
wiggler  ultimately  breaks  down  in  this  limit. 
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Fig.  2  —  Schematic  representation  of  the  dispersion  curves 
the  cyclotron  and  the  pure  electron'll 
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Fig.  3  —  Graph  of  <t>  versus  j80  for  both  group  I  and  group  il  orbiis  for  a  beam  energy  of  1.5  MeV  and  a  wiggler 

amplitude  and  period  such  that  fl  „/*„.<•  —  05 


21 


Fig.  4  (Continued)  —  Schematic  representation  of  the  interaction  for  group  I  (a)  and  group  II  < b)  orbits  when  d>  >  0 

The  dotted  line  denotes  complex  conjugate  roots. 
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Re  k  +  /kw 

Fig.  10  —  Graph  of  versus  Re*.  *„  lor  group  II  ('!>  <  0)  orbits  and  “  I 
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